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Abstract 

We re-examine the Boussinesq hypothesis of an effective turbulent viscosity within 
the context of simple closure considerations for models of strong magnetohydrody- 
namic turbulence. Reynolds-stress and turbulent Maxwell-stress closure models will 
necessarily introduce a suite of transport coefficients, all of which are to some degree 
model-dependent. One of the most important coefficients is the relaxation time for 
the turbulent Maxwell stress, which until recently has been relatively ignored. We 
discuss this relaxation within the context of magnetohydrodynamic turbulence in 
steady high Reynolds-number, high magnetic-Reynolds-number shearing flows. The 
relaxation time for the turbulent Maxwell stress is not limited by the shear time 
scale, in contrast with Reynolds-stress closure models for purely hydro dynamical 
turbulence. The anisotropy of the turbulent stress tensor for magnetohydrodynamic 
turbulence, even for the case of zero mean-field considered here, can therefore not 
be neglected. This shear-generated anisotropy can be interpreted as being due to 
an effective turbulent elasticity, in analogy to the Boussinesq turbulent viscosity. 
We claim that this turbulent elasticity should be important for any astrophysi- 
cal problem in which the turbulent stress in quasi-steady shear has been treated 
phenomenologically with an effective viscosity. 
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1 Introduction 



The Boussinesq hypothesis of an effective turbulent viscosity is still widely 
used both in engineering appUcations and in astrophysics. This is true despite 
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the passage of over one hundred years since the birth of the concept. There 
are a very large number of models for turbulent stress that give substantially 
better results than a simple effective viscosity, and in some situations, such as 
strongly time-dependent flows, the use of an effective viscosity gives notori- 
ously bad results. Despite this, the concept of an effective turbulent viscosity 
has proven to be so persistent because it is useful, simple, and because it 
is based on an appealing analogy between the random motions of turbulent 
blobs and the kinetic theory of thermal motions of molecules. Thus, effec- 
tive turbulent viscosity is commonly invoked as a good crude "lowest-order" 
approximation, for problems in which the Reynolds number is very large. 



For problems in which the magnetic Reynolds number is also very large (as is 
often the case in astrophysics), we ask if a picture of turbulence as colliding 
blobs of fluid is at all useful even as a crude approximation. A more appropri- 
ate analogy for the turbulent component of the field may be not the kinetic 
theory of point-like molecules, but rather the kinetic theory of long polymers 
in melt or in solution ( |Ogilvie| , |2001| [William^ , |2001| ). For the crudest possible 
model, such an analogy would suggest the concept of an effective turbulent 
elasticity, in addition to the effective turbulent viscosity. Thus, corresponding 
to the Boussinesq approximation in hydrodynamic turbulence, the correspond- 
ing "lowest-order" approximation for the stress in MHD turbulence should in- 
clude an additional coefficient, in analogy to a simple viscoelastic fluid. Just as 
the inclusion of elasticity in the stress response of an ordinary laboratory fluid 
can have dramatic consequences, we suggest that any astrophysical problem 
in which turbulence has been treated phenomenologically with the Boussinesq 
viscosity might have dramatically different solutions if an effective elasticity 
were included as well. 



In particular, for the purposes of quasi-steady shear, we claim that the anisotropy 
of the stress tensor is of much greater significance in the case of MHD turbu- 
lence than in the case of purely hydrodynamic turbulence. Thus, while inves- 
tigations of the idealized problem of homogeneous istotropic hydrodynamic 
turbulence may lead to a greater understanding of the dynamical importance 
of turbulence to the large-scale mean flow of a turbulent fluid, we expect that, 
quite broadly speaking, models for MHD turbulence that are to be of some 
use in predicting the behavior of the mean flow must address directly the 
anisotropy of the turbulent Maxwell stresses, irrespective of the presence or 
absence of any mean field. 



Many of the more recent simulations of MHD turbulence in a shearing en- 
voronment (our main interest in this paper) are simulations (e.g. [Hawley et al 



(|1995|) ) of the Velikov-Chandrasekhar-Balbus-Hawley magnet orot at ional in- 
stability (MRI), and we will refer to the results of these simulations for guid- 
ance. See also |Williams| ( |2002|) for a comparison of simulations of the MRI 
with some viscoelastic models. 



2 



For our purposes, the dominant feature of simulations of the MRI is that the 
largest component of the turbulent stress tensors is not the "viscous" r6 com- 
ponent, but rather the streamwise 66 component. (Note that in simulations of 
the MRI, 6 is the direction of the shear, and r is the cross-shear direction.) In 
fact, our original reason for investigating viscoelastic models was as a mecha- 
nism to coUimate and possibly drive jets using this 66 component of the stress 
tensor. 

Most importantly we wish to point out that the interpretation of the turbulent 
stress as an effective viscous stress can lead to qualitatively erroneous predic- 
tions. In view of the fact pointed out above that the "viscous" component is 
in fact not even the largest component of the stress tensor, a prescription for 
the turbulent stress tensor n*"^^ of the form 



where /xturb and /iturb effective shear and bulk viscosities respectively, is 
untenable, even as a gross approximation. 

One possible approach to the problem is to generalize the viscosity to a tensor 
viscosity. This approach ignores the fact, however, that the large azimuthal 
component of stress in MRI-driven shear turbulence is fundamentally due to 
the vector advection properties of the induction equation. 



2 Notation 

As is conventionally done in the Reynolds decomposition in hydrodynamic 
turbulence, it is convenient here to decompose all fluid variables into mean 
and fluctuating, or "turbulent," parts. The former will be denoted with an 
overbar, and the latter with a prime. Thus, for example, Bi = Bi + B[. We 
employ standard ensemble rules of averaging, so that averaging commutes with 
time and spatial derivatives. It is also useful to think of this averaging as, say, 
a time or spatial averaging over some convenient timescale or lengthscale, 
although strictly speaking this is incorrect, because such averaging introduces 
additional terms. 

Furthermore, for the sake of definiteness, we will make liberal use of Carte- 
sian index notation for vectors and tensors, although we may occasionally 
use abstract vector notation when there is no chance of ambiguity. To make 
lengthy equations a bit more palatable, we denote symmetrization and an- 
tisymmetrization on indices with curly and square brackets respectively, so 



turb 
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that, e.g., 

A{ij} = Aij + Aji. 



The generahzation of the given equations to curvihnear ( "covariant and con- 
travariant") form is straightforward. 

We use Heaviside-Lorentz units for the magnetic field. The use of Gaussian 
units and the attendant factors of 1/ v^ivr would make lengthy equations even 
longer, and this would detract from readability. The conversion to Gaussian 
units is 

j^[Gauss] /~, j^[Heaviside—Lorentz] 



The full Maxwell stress tensor, under the approximation that \E\ ^ \B\, has 
components BiBj — ^B'^6ij. The turbulent Maxwell stress is the Maxwell stress 
due to the turbulent component of the field, B'. We find it more convenient to 
work with the magnetic cross-correlation tensor Mij where Mjj = B[B'j. Then 
the turbulent Maxwell stress is 

M,. = M,. - MkkSij. 



3 Invariant Tensor Advection Operators 

In perfect flux-freezing the stress tensor Mij = B'^Bj evolves by the action of 

the (modified) upper-convected invariant derivatve such that D\^\M)ij = 
0, where 

(Dl"lM)i,- = {dt + Vkdk)Mi, - {dkiH)Mkj - M,k{dkVj) + 2{dkVk)M,j. (2) 

(The required modification of the conventional upper-convected invariant deriva- 
tive is the inclusion of the term proportional to the divergence of the velocity). 
This relationship is simply a direct result of the vector advection equation for 
the magnetic field in perfect MHD, and does not rely on breaking the field 
into mean and fluctuating parts, although we have done so here. This equation 
shows how, for example, the stress Mij is distorted by shear (see figure 1). 

Generally speaking, tensors formed from a dyad of vectors that obey vector 
advection equations will often produce a tensor advection equation that is 
similar to the above: 
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Consider the scalar advection operator 



dt + Vkdk- 



This may be generahzed to a vector advection operator in more than one way. 
Note that the vector advection operator acting on either the field B in MHD or 
the vorticity uj in simple hydrodynamics is, in the incompressible case (writing 
the vector as w), 

dtWi + VkdkWi - WkdkVi, 

refiecting that these vectors are stretched and aligned with the shear. In con- 
trast, in the Reynolds decomposition, the advection operator acting on the 
turbulent velocity field v'^ is 

dtv'i + Vkdkv'i + v'^dkVi + v'l^dkVi 

The last term is the origin of the famous triple-correlation; here we focus on 
the penultimate term, which is similar to the last term in expression (3) except 
for a change of sign. Certainly more complicated examples are possible. For 
instance, the advection of the gradient Qi = dicj) of a passive scalar is 

dtQi + VkdkQi + QkdkVi - eijkQjUJk = 

In all cases above, however, there is a term WkdkVi that, depending upon 
whether it appears with a plus sign or a minus sign, tends to anti-align or 
align the vector Wi with the shear. Thus, if the stress tensor is formed from 
vectors that obey such equations, the stress will necessarily be distorted by the 
Lagrangian deformation of the fiuid. Upper- and lower-convected derivatives 
are two tensor generalizations of the advection operator that arise when the 
tensors are formed from vectors that are aligned or anti-aligned by the shear. 

One finds then that, for example, the Reynolds stress is distorted by the 
background shear by the advection operator (the lower-convected invariant 
derivative) which, in the incompressible case, may be written as 



The second and third terms in this operator definition have long been rec- 
ognized as being quite important to the theory of hydro dynamical turbulent 
transport processes. These terms distort the Reynolds stress, and are as well 
a source for the Reynolds stress. The distortion that shear induces on Rjj 
is opposite to that induced on Mij, and in such a manner as to produce a 




(3) 
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Fig. 1. The turbulent magnetic field lines are distorted by shear. In equilibrium, 
the degree of anisotropy is determined by the balance of production, distortion, and 
dissipation, as described below. 

"viscous" -like component to the stress. As the turbulent Maxwell stress and 
the Reynolds stress appear with different signs in the turbulent stress tensor, 

n7^ = -R,,+M,„ (4) 

the signs of the viscous component of the Maxwell stress and the Reynolds 
stress are the same. It is of course just such a distortion of the stress tensor by 
the deformation of the fluid (including shear) that is the origin of molecular 
viscosity, where in that case the stress tensor is formed from the sums of the 
contributions mViVj to the stress tensor from individual particles. 



4 Reynolds Decomposition 

It is helpful to write down in full the Reynolds decomposition for MHD tur- 
bulence. Here let us concentrate on the const ant- density case, p = p. The 
continuity equation and "magnetic Gauss's law" are then 

dkVk = dkv'k = 
dkB„ = d^B'^ = 

The mean-field velocity transport equation is 

pdtVi + pVkdkVi = -diP + BkdkBi - \di{B'^) - 
dkRik + dkMik - \diMkk + p^dkkVi 
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The unmodeled transport equation for the Reynolds stress is obtained after 
some tedious but straightforward algebra: 



dtRij + VkdkRij + {dkVi)Rkj + Rik{dkVj) + pdk{v[v'jV'f,) = 

-4 {d^yiBkB',)) - lv[^ {d^yiB'.B',)) + /i {dlR., - 2idkvi)idkv'^)) (5) 

There are three magnetic curvature terms, namely those arising from B' ■ V-B, 
those from B ■ VB', and those from B' ■ V-B'. Likewise, there are two magnetic 
pressure terms, namely those arising from the cross-correlation pressure BB' 
and those from the pure turbulent field pressure B'B'. It should be noted 
that in all terms in which both the velocity and the field appear, the velocity 
appears once and the field appears twice. 

We write the induction equation as 

dtBi + VkdkBi - BkdkVi + Bi{dkVk) = vdlBi. 

From the induction equation we find that the mean field evolves according to 
dtBi + VkdkBi - BkdkVi = Sfc^J^^ + r]dlBi 

where we have defined the antisymmetric matrix 
(It is also convenient to define the symmetric matrix 

although we do not use this here.) From this we can see that the only sources 
for a mean field B are a stretching component (which also appears in the 
induction equation itself) and a turbulent source term which is more often 
written as 

dkA^^ = X (V^)) ^ 

As we are mainly interested in the zero-mean-field case, we do not discuss the 
mean-field induction equation any further. 

The turbulent Maxwell-like stress Mij obeys the transport equation 
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+ry {dlM,, - 2{d,Bi){d,B'^)) (6) 



We have so far discussed the Reynolds stress and the turbulent Maxwell stress, 
but a word is in order regarding other approaches to turbulence. In particular 
we wish to address briefly the dynamics of turbulent vorticity. It is tempting 
to explore the apparent symmetry in the transport of the vorticity uj and 
the magenetic field B. In high Reynolds-number flow the vorticity tends to 
concentrate on discrete vortex lines. In the extreme case of a superfluid the 
vorticity of the fluid is identically zero except on discrete quantized vortex 
lines, just as in a superconductor the magnetic field is zero except on discrete 
quantized field lines; there is certainly an interesting symmetry to be expolored 
in the interaction of discrete vortex lines and discrete field lines in the case 
of turbulence in a superconducting superfluid. One may suspect that this 
analogy will also hold in the case of high Reynolds-number, high magnetic 
Reynolds-number turbulence. 



However, this apparent symmetry is far from perfect. Neither vector is trans- 
ported as a passive field. The magnetic field exerts its influence on the bulk 
flow through the appearance of the divergence of the Maxwell stress as a bulk 
force in the fluid equations. In contrast, the vorticity does not so much in- 
fluence the flow as it does determine it, insofar as the velocity is determined 
directly from the vorticity through the Biot-Savart relation, and conversely 
the vorticity is simply the curl of the velocity. As a result, a vortex line, if 
bent, can "wrap itself up" by the self-action of its own vorticity, unlike a field 
line. 



So, to the extent that vorticity dynamics may be used - in analogy to polymer 
dynamics in solution ( Phorin| , [L988|) - as a basis for turbulence models and 
quasi-viscoelastic behavior in particular, we expect that the corresponding 
analogy of turbulent field lines to polymers should be much more fruitful. Note 
also that a stress tensor may be formed directly from the field, whereas the 
relationship between the turbulent vorticity and the Reynolds stress is not 
straightforward. For this reason we do not treat the dynamics of turbulent 
vorticity and we resort entirely to simple Reynolds stress modeling for the 
stress due to turbulent hydrodynamic motions. 
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5 Boussinesq-Like Approximations 



5.1 Simplest Case: No coupling 



An ordinary differential model, such as the viscoelastic Maxwell model, con- 
tains in principle all that is needed in order to determine the stress in a ffuid, 
given initial conditions and a rate-of-strain history. In steady, so-called "curvi- 
lineal" ffows with closed or nearly-closed streamlines (such as the azimuthal 
Stokes shear ffow in the neighborhood of a rotating sphere), it is sometimes 
useful to rewrite the solutions of these differential equations in closed form. 

For example, the upper-convected Maxwell model may be written in closed 
form as 

t 

^ = Y2 J (Cr'(r)-l)exp(-(t-r)/A) dr, 

— oo 



where ^ is the Finger tensor (see [Joseph] (|1990| ) p. 14). Regardless of whether 



an integral expression may be found for a given model, another closed-form 
way of writing the stress for steady, curvilineal flows is to expand the stress as 
a series solution in the powers of the (symmetrized) spatial derivatives of the 
velocity. The conventional way of doing this is to make use of powers of the 
Rivlin-Eriksen tensors A'"', defined below. The first Rivlin-Eriksen tensor, A^^\ 
is the symmetrized velocity gradient (defined in terms of the mean velocity Vi), 
so, with the addition of a "pressure" term proportional to 6ij, the Boussinesq 
model may be written as 



for incompressible flow. (If the flow is compressible, it is conventional to break 
the rate-of-strain tensor into the divergence and a traceless part as we have 
done ealier; the coefficient of the former is then the turbulent bulk viscosity.) 

A Boussinesq-like approximation may be "derived," in a very crude sense, from 
simple turbulence models A full turbulence model should make considera- 
tion of the source of the turbulence, such as an instability such as convection 
or the MRI. However, we point out that a simple turbulent viscosity works 
well as a first approximation when solving for the turbulent transport of flow 
through the lumen of a pipe (i.e. pipe Poiseuille flow), for which there is no 
known instability, linear or otherwise. 

^ The author is grateful to Michael Steinkamp for pointing this out to him. 
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The unmodeled Reynolds transport equation for incompressible fluid, under 
conventional rules of ensemble averaging, and ignoring effects such as rotation 
and thermal convective driving, is 

dtRij + VkdkKij + {dkVj)Rik + {dkVi)Kkj = 
-dkWH) - + W^) + ^ [^l^^^ - '^{dkv[){dkv'^)) . (7) 

(This is simply eqn. (5) rewritten without the terms in which the magnetic 
fleld occurs.) The left-hand side of this equation represents the advection and 
distortion of the Reynolds stress by the mean flow. The right-hand side of 
the equation is the hard part. The crudest models for the right-hand side 
of this transport equation will include terms proportional to 5ij and terms 
proportional to Rjj. Effectively, we have advection and distortion on the left 
hand side (this distortion is sometimes referred to as a source, since it changes 
the total kinetic energy iRfc^), and on the right hand side we have a source 
proportional to 5ij and a sink proportional to R.^. It is somewhat abusive to 
call these terms source and sink terms, but the point is that the one typically 
has a positive coefficient and the other typically has a negative coefficient. For 
example, consider the Rotta return-to-isotropy term, which is a model for one 
of the effects of the correlations of turbulence with pressure fluctuations, i.e. 
the second of the three terms on the right-hand side of the above equation 
(note that pressure fluctuations also carry sound away, but this is not included 
in the Rotta term). This may be written as, e.g., 

where e and K are the turbulence dissipation rate and the turbulent kinetic 
energy respectively, and the Rotta coefficient Cr is a positive constant 0. 
Likewise the viscous part of this transport equation — i.e. the last term on 
the right-hand-side — has a diffusive part, which is usually ignored, and a 
dissipative part that is most simply modeled as being proportional to the 
Reynolds stress R^j itself. 

We thus write 

dtRij + VkdkRij + {dkVj)'R.ik + {dkVi)'Rkj = a^Sij - s'^Rij 

where and are numerical coefficients (which may depend upon certain 
characteristics of the flow, such as the shear, the buoyant convective driving, 
and so forth). Note that is a relaxation time. We may then solve for Rij, in 

^ Francis Harlow, personal communication 
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the case of steady shear (so that the advective operator dt + Vkdk gives zero), 
by plugging the above equation back into itself to obtain an infinite series, 

Ti - ^ n A..-«2^ /lW_L.3^ ^[3] I 

This may be seen from the following relation for A'"': 

In fact, it is not at all necessary to assume that the "source" term is propor- 
tional to If we instead write 

dtRij + VkdkRij + {dkVj)'Kik + {dkVi)'Rkj = slf^ - s~'^Kij 

where Slf^ is some unspecified tensor, we obtain 

R,- = s^sff - slUS^%, + slL^{S^%, - sll.'{S^%, + ... 

where L is the linear operator 

= {dkVijSkj + S^^\dkVj), 

so we may formally write 

If the relaxation time is much shorter than the shear time scale (which we 
denote l/7s), then one expects the series to converge rapidly. 

Despite our formal solution above, let us concentrate on the original case, 
S'lj^^ = a^Sij. To first order in the shear one obtains then 

^ij = ^R^^R^ij - ^R^^R^diVj + djVi) = Pturb - ^turh{diVj + djVi). 

The exact form of the Reynolds stress in this case, where s^7s <C 1, is actu- 
ally quite sensitive to the source term, which we have here modeled as being 
proportional to Sij. Nevertheless, the approximation that one obtains — an 
effective turbulent pressure and an effective turbulent viscosity — has proven 
quite useful over the past century. 

The product s^'-fs, which one might call a "turbulent hydrodynamic Weissenburg- 
Deborah number," is in fact generally less than or of the order of one in pure 
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hydrodynamic turbulence. One can make the argument in very rough fash- 
ion using eddy-viscosity type arguments. Colhding turbulent blobs lose their 
identity on the order of a coherence timescale, which is on the order of an 
"eddy turnover" timescale, which is of the order of the turbulent lengthscale 
^turb divided by the turbulent velocity scale Wturb- The magnitude of the differ- 
ence in the corresponding background mean flow from one side of the eddy to 
the other is ituAls- If this velocity is greater than the turbulent velocity, the 
baackground shear will rapidly disrupt the eddy. Thus iturhls < ^^'turb- But this 
implies that 5^7^ < 1, if ~ •^turb/'^^turb- Likewise, in the case of the contin- 
uum approximation in hydro dyanmics, the molecular collision frequency must 
exceed the shear rate for the concept of a Stokes molecular viscosity to be 
terribly useful. 

The admittedly very crude argument given above does not hold at all in the 
case of a high magnetic-Reynolds-number turbulent magnetic field, however. 
In perfect fiux-freezing, a magnetic field line maintains its identity, or "coher- 
ence," indefinitely. 

In reality, dissipative effects will become important at some point. This is 
true even in ideal MHD, we expect, if we extend the word "dissipative" to 
include the effects of turbulent cascades and inverse cascades, which will take 
energy on intemediate scales and transport that energy either to very small 
scales or to larger scales. (Whether energy that is transported to small scales 
is ultimately dissipated in the molecular sense is not necessarily relevant.) It is 
not at all clear, however, that these quasi-dissipative effects become important 
on a shear time scale. 

Consider the transport equation for the turbulent magnetic cross-correlation 
tensor Mij as written in eqn. (6). If we write 

we obtain, in a similar matter to the above. 

Again, taking the simplest case, S\j'^^ = a^^jSij, we obtain 
Mij = s,,ajij + sl^a^^{diVj + djVi) + s^^a^^A^j + ... 

Note that the signs of the coefficients here do not alternate, so that if the 
relaxation times were equal, = then the series solution for M might 
not be expected to converge as rapidly as that for R. In reality, however, we 
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note that for simple shear flow (or "curvlineal" flow), both series converge 
exactly with only three terms, as A^"! = for n > 2. 

The full turbulent stress tensor is 

Hij = — Rij + Mij — ^Mi:k6ij, 

recalling that the Maxwell stress tensor 7-[A^«^"'eZ/] ^^-^ |^| ^ |^|^ 

[Maxwell] _ r p 1 



SO that the turbulent Maxwell stress tensor Mj^ is 



We may then write, again for the simplest case where S)^^ = a^6ij and Slf^^ 
djij^ij. that 



The identification of the term proportional to 6ij above as a turbulent pressure 
is not unique, as the rest of the terms are not all traceless. This ambiguity nat- 
urally arises from the identification of as the magnetic pressure. Consider 
the case M^j = mSij. Then the turbulent Maxwell stress is 

3 

Mij = Mij - \Mkk5ij = m6ij - -m5ij = -\m5ij 



so that the Maxwell stress is proportional to 6ij, which looks like a pressure, 
except that the nominal magnetic pressure contribution to the stress is actually 
— (3/2)m5jj. With this caveat in mind, then, we identify the term proportional 
to 5ij as a pressure. This enables us to write 

^ij = ~Ptuvh^ij + /^turb^y + Cturb^lj ) (8) 



where Cturb is the effective turbulent elasticity. (One alternative choice for 
defining the turbulent pressure is to assume that the pressure is proportional 
to the trace of the full stress tensor, so that 

n - P A 4_ „ 4W-U^ /^4[2] -^4[2]a V 

J-J-ji — turbOij ~r /^turb^jj "T Sturb \ ~ '^^kk^ij ) ' 
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it is not immediately clear if this is a preferable definition or not.) We may 
then define an effective turbulent relaxation rate as 

Cturb 

Seflf = 

/^turb 

and an effective turbulent Weissenburg number 

Weturb = IsSeS 

The elastic term might also be thought of as a higher-order viscous term, for 
this restricted case of steady curvilineal flow. 

5.2 Inclusion of coupling between Reynolds and turbulent Maxwell stresses 

In fact, we also obtain an expression of the form (8) if we take the more general 
case that 

dtRij + VkdkRij + {dkVj)Rik + {dkV,i)Rkj = ajij + h^uMij - s^^Rjj 
dtMi^ + VkdkM.j - {dkVi)Mkj - {dkVj)Mik = Sjf ' + ^mrR^,- - (l/s.J^u 

That is to say, in this case, we make use of the available tensors 6ij, Rjj, 
and Mij. The new terms (with coefficients 6^^^ and b^^^^) represent exchange 
terms. As such, the coefficients must have opposite signs. Such a form may be 
motivated on the general grounds that energy should be exchanged between 
the Reynolds stress and the turbulent Maxwell stress; further motivation is 
found through a more careful analysis of the transport equations as discussed 
below. 

Let us rewrite the Mij transport equation, highlighting a few terms, and re- 
placing the resistive term with a simple dissipative term. We have 

dtMij + VkdkMij - Mik{dkVj) - {dkVi)Mkj+ = 
-{dkB^B^, + BkW^[jB^^ - duiMBWk) + (^^^4PPI - (9) 

Terms written in light grey are those in which the mean field Bi appears. The 
terms in green exchange energy between Rjj and M^j through the respective 
terms v'^dkB[ and B^OkV^ that appear in the induction equation. (It might 
seem that we should discuss energy exchange between the Reynolds stress R 
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and the Maxwell stress M rather than between R and M. However, note that 
the turbulent energy density f/turb is 

t^urb = ^R-fc/c — |Mfcfc = |Rfcfc + ^Mkk, (10) 



so that we recover that the turbulent field energy density is (5^5^) /2, as it 
should be.) 

Similarly, for the Reynolds stress, we write 



dtRij + VkdkRij + {dkVi)Rkj + Rik{dkVj) + pdk{v[v'jv'k) = 

-v[^ {d^^iBkBi)) - lv[^ {d^^iBiBi)) - s-'R,, (11) 
Again, R-M exchange terms are highlighted in green. 

One would expect the exchange terms to be proportional to the mean of 
v' ■ {B' ■ VB' — iV(-B')^); this is easily confirmed. The trace of the (zero mean- 
field) exchange terms in the Reynolds stress transport equation is 

Tr{terms) = 2v',B',d,B', - v',d,{B'^B'^) (12) 

whereas the trace of the exchange terms in the Maxwell-like (i.e. M) stress 
transport equation is 

Tr(terms) = 2v[B',d,B', - v[d,{B',B',) + 2du 

this last term being in conservative form and so identically zero in the homo- 
geneous case. 

In a more sophisticated model, the exchange between Rjj and Mij might 
be modeled through a transport equation for the triple-correlation exchange 
terms described above. Here, however, we simply model this exchange in the 
transport equations for Rjj and Mij by assuming that the exchange is propor- 
tional to these tensors themselves. 

With this assumption, R^j and Mij have solutions of the form 

Rii = ro5,j + riAjJl + rs^gl 
Mij = mo6ij + miAf^ + m2Afj 



(13) 
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So, the expression (8) is again obtained. Specifically, the solution is 

^^0(1 ^ ^R^RM^M^Mr) ~ ^M^M + ^M^MR^R^R 

5.3 Effects of Rotation and Stratification 



It might be expected that rotation is not something to be put into a turbu- 
lence transport model; it ought to exist already in the trasport model, so that 
rotation is just a matter of doing the proper coordinate transformation. This 
philosophy is also in keeping with the principle of material frame indifference, 
according to which a fluid should not know (at least as far as its constitutive 
relations are concerned) about the overall rotation of the system. Actually 
the issue is a bit more subtle than this, since the form of the Reynolds and 
Maxwell stress transport equations described above (or, typically, in their un- 
modeled state) will be invariant with respect to transformation to a rotating 
reference frame, but this means that the action of the Coriolis force on the 
turbulent blobs in the rotating reference frame is ignored (see, e.g., the dis- 
cussions in Pope| (|2000|) and in pR.iidiger| ( |1989|) ). Inclusion of the Coriolis force 



may appear as the introduction of new terms with symmetrized products of 
the rotation tensor VLij = eijk^k, as described below. Such considerations are 
important for a model of the MRI, as well as for the more simple case of 
Taylor instability when the Rayleigh stability criterion is not satisfied. The 
inclusion of the effects of rotation is an ongoing project that we intend for a 
future paper. These considerations must be addressed in order to reproduce 
the result from simulations that the rr component of the stress is larger than 
the zz component of the stress, not to mention the presence of the MRI itself. 

It would also be nice to address the effects of stratification on the turbu- 
lence, such as is a nessesary ingredient for the treatment of convective MHD 
turbulence. Consideration of turbulent elasticity should be important for the 
turbulence in the solar convective zone (SCZ), for example, even though the 
turbulence in the SCZ is not driven (primarily, at least) by MHD instabili- 
ties but rather purely hydrodynamical ones. For incompressible but variable- 
density stratified fluids in the laboratory, one finds an additional driving term 
in the Reynolds stress transport equation that is proportional to the product 
of the pressure gradient and the density gradient. In cases of astrophysical 
interest, where compressiblity becomes important, the important considera- 
tion is the direction of the entropy gradient (or, the gradient in the potential 
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temperature). This will necessarily introduce anisotropics into the turbulent 
stress tensor, just as a consideration of the sources of the MRI will produce 
anisotropics that are not adequately addressed by the considerations outlined 
in this paper, as mentioned above. Inclusion of these effects may appear as the 
introduction of new terms with symmetrized products of the vectors gi and 
diP and dip, or other thermodynamic derivatives, depending upon the choice 
made by the modeler for the independent thermodynamic quantities. Again, 
a proper inclusion of these effects must follow from a model that is built from 
the transport equations in which — unlike what we have presented here — 
the assumption of incompressibility is relaxed. This is also an ongoing project 
that we intend for a future paper. 

We can, however, comment briefly upon the combined effects of rotation and 
stratification as they manifest themselves in stellar convective zones through 
the so-caled "A-effect" ([Rudiger| , |1989| ). This introduces additional terms that, 
in the SCZ, are to lowest order proportional to the dyad product {g x Vl)g. 
(Note that inclusion of the rotation VLi violates one of the assumptions that 
forms the basis of the expansion in terms of Rivlin-Eriksen tensors, namely 
material frame indifference, but again, this is to be expected for a turbulence 
model that takes account of Coriolis forces, as mentioned above.) 

For the A-effect as described by [Riidigei] ( |1989| ) (in particular, section 4.6), 
the additional terms appear in the prescription for the equilibrium Reynolds 
stress (as opposed to the Reynolds stress transport equation), and for the 
lowest-order A-effect the symmetrized product mentioned above, 

j^lambda _ -^^ {g^gk^lkj + ^ikQugj) + • • • , 



produces additional components Rj.^" and, to higher order, R^*^ . These 
decouple from the elastic terms in the stress, so that the total turbulent stress 
tensor may be written as 

TT pturb I TTvisc I TTelastic _i_ Trlambda 

J-J-ij — -T jj "T 11 jj "I- lljj "I- lljj 



The form of the turbulent stress tensor may then be written as an expression 
of the form (compare [Rekowski and Riidigei) ( p.998| ) eqn. 17) 

= -P5ij + pNijke{dkVi) + pZijkemn{dkVi){d,nVn) - pKjk^k (14) 



so that the various effects are simply additive. 
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6 Conclusions 



We have discussed some modifications to the basic Boussinesq ansatz that 
are motivated by an analogy between the turbulent component of the mag- 
netic field in MHD turbulence and the dynamics of polymers in solution. This 
is part of an ongoing project begun by us to study MHD turbulence begin- 
ning with the Reynolds decomposition of the field, with particular focus on the 
small-scale tangled component in preference to the large-scale mean-field com- 
ponent. It is our belief that undue focus on the latter, apparently for the sake 
of understanding dynamo behavior, has obscured the important dynamical 
effects that the former may have in practical MHD turbulent flow problems, 
and that furthermore the Boussinesq effective viscosity had been adopted for 
such flows without sufficient consideration as to whether the initial motiva- 
tion for an effective turbulent viscosity carried over from hydrodynamics to 
MHD. While a viscoelastic picture of MHD turbulence, built from — or at 
least inspired by — a polymer analogy, may not be "correct" in some absolute 
sense (as pointed out by a previous anonymous referee), we feel that it is at 
least useful, which is all that a model can ask to be. Of course, there are many 
effects that naturally arise from viscoelastic models, as described at length in 
any textbook on the matter, which we hope to discuss in future work. For 
example, the turbulent medium supports Alfven waves that can be treated 
easily within a viscoelastic picture; these turbulent Alfven waves have a com- 
plex dispersion relation, reflecting the dissipative nature of turbulence. These 
waves do not appear at all in a purely viscous model. We take this as a further 
indication that a viscoelastic picture of MHD turbulence may be useful. 

Here we have focused on one aspect common to any viscoelastic model of MHD 
turbulence, which is to say any model that takes into account the distortion of 
the turbulent Maxwell stress by the inclusion of the proper tensor derivatives 
and a relaxation time, namely the appearance of signiflcant anisotropy of the 
turbulent stress tensor in the direction of the shear. This effect should be 
important in any shearing environment in MHD turbulence, including but not 
limited to the MRl (and in particular the inner regions of accretion disks and 
the jet-launching process) and to the turbulence in stellar convection zones. 
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